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ABSTRACT 

Previous analysis of the fossil-group/cluster RXJ1159+5531 with X-ray observations front a central 
Chandra pointing and an offset-North Suzaku pointing indicate a radial intracluster medium (ICM) 
entropy profile at the virial radius (i? v ir) consistent with predictions from gravity-only cosmologi¬ 
cal simulations, in contrast to other cool-core clusters. To examine the generality of these results, we 
present three new Suzaku observations that, in conjunction with the North pointing, provide complete 
azimuthal coverage out to R v ; r . With two new Chandra ACIS-I observations overlapping the North 
Suzaku pointing, we have resolved >50% of the cosmic X-ray background there. We present radial 
profiles of the ICM density, temperature, entropy, and pressure obtained for each of the four direc¬ 
tions. We measure only modest azimuthal scatter in the ICM properties at R 200 between the Suzaku 
pointings: 7.6% in temperature and 8.6% in density, while the systematic errors can be significant. 
The temperature scatter, in particular, is lower than that studied at R 200 for a small number of other 
clusters observed with Suzaku. These azimuthal measurements verify that RXJ1159+5531 is a regu¬ 
lar, highly relaxed system. The well-behaved entropy profiles we have measured for RXJ1159+5531 
disfavor the weakening of the accretion shock as an explanation of the entropy flattening found in 
other cool-core clusters but is consistent with other explanations such as gas clumping, electron-ion 
non-equilibrium, non-thermal pressure support, and cosmic ray acceleration. Finally, we mention that 
the large-scale galaxy density distribution of RXJ1159+5531 seems to have little impact on its gas 
properties near R V1I . 

Subject headings: X-rays: galaxies: luminosity - galaxies: ISM - galaxies: elliptical and lenticular 
Clusters of galaxies: intracluster medium 


1. INTRODUCTION 

Galaxy clusters represent the final stage of the hier¬ 
archical formation, and they are potent laboratories for 
testing models of structure formation. They are espe¬ 
cially valuable when considering global quantities com¬ 
puted within their virial radii (R v i r ), since the outskirts 
of galaxy clusters contain most of their baryons, dark 
matter, and metal content. The entropy profile is a sen¬ 
sitive indicator of non-gravitational heating of the intra¬ 
cluster medium (ICM), since entropy is conserved in an 
adiabatic process. In studies of galaxy clusters, entropy 
is conventionally defined as K(r) = T(r)/n(r) 2 / 3 where 
T and n are gas temperature and density. The deviation 
of the entropy profile from that expected from pure grav¬ 
itational collapse K{r) cx r 11 (Voit et al. 2005) reflects 
the role of non-gravitational processes in cluster forma¬ 
tion. Cluster outskirt near the virial radius are the front 
lines of cluster formation and are still growing. Primor¬ 
dial gas and sub-structures have been continuously ac¬ 
creted into the ICM along large filamentary structures, 
leaving a shock-heated region near the R v ; r . 

To date, a dozen of massive galaxy clusters have been 
observed with Suzaku out to 7? v ir- These studies found 
that the gas properties near I? v j r disagree with the predic¬ 
tions of gravity-only cosmological simulations (e.g. Su et 
al. 2013, Bautz et ah 2009; George et al. 2009; Simionescu 
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et al. 2011; Hoshino et al. 2010; Walker et al. 2012a; 
Walker et al. 2013, see Reiprich et al. 2013 for a review). 
In particular, outside of -Rsocj^the observed entropies are 
significantly less than the gravity-only prediction. This 
entropy deficit cannot be explained by feedback, the ef¬ 
fect of which is expected to be small at such large radii, 
and should instead produce an excess of entropy over 
that produced by only gravitational evolution. 

In the case of the Perseus Cluster, Simionescu et al. 
(2011) also measured an enclosed gas fraction within 
I ?200 that exceeds the cosmic value by 50%. Simionescu 
et al. (2011), in particular, advocate a clumpy ICM to 
reconcile the observations with structure formation mod¬ 
els, since clumped gas emit more efficiently than uni¬ 
formly distributed gas. The clumping factor, defined as 


with C > f, is used to describe the deviation from uni¬ 
formly distributed gas. In the case of the Perseus Cluster, 

2 Ra is the radius within which the average density is A times 
the cosmological critical density. R 200 ~ Rv ir am! R 500 ~ 0.6 R vlr . 
We adopt -R 200 as i? v ir the fiducial virial radius in our study for 
two reasons. First, R 200 is most commonly used as in X-ray clus¬ 
ter studies. Second, we are able to obtain interesting constraints 
on the ICM properties out to that radius with all four Suzaku ob¬ 
servations. Since .Rios was previously employed by Humphrey et 
al. (2012) to study the ICM properties of RXJ1159+5531, we also 
make some comparisons of the ICM properties at that radius. 
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a clumping factor of ~16 is required at i? 200 - C^7 and 
~9 are required for Abell 1835 (Bonamente et al. 2013) 
and PKS 0745-191 (Walker et al. 2012b) respectively to 
reconcile the observed entropy profile and the expected 
r 11 power-law model from gravity-only simulations. In 
contrast, significantly smaller clumping factors are pre¬ 
dicted for simulated clusters. For example, Nagai & Lau 
(2011) obtained C =1.3-2 at the R V1I . Similarly, Vazza 
et al. (2013) and Zhuravleva et al. (2013) report C <3 in 
the cluster outskirts. 

Apparently, other factors besides gas clumping need 
to be taken into account to explain the observed en¬ 
tropy profiles. Hoshino et al. (2010) and Akamatsu et al. 
(2011) attribute the entropy flattening to electrons and 
ions being out of thermal equilibrium due to recent accre¬ 
tion shocks (Hoshino et al. 2010, Akamatsu et al. 2011); 
however, some simulations suggest such effects should 
not be significant (Wong & Sarazin 2009). Femiano & 
Lapi (2014) proposed that the rapid radial decrement 
of the temperature caused by non-gravitational effects is 
responsible for the observed entropy flattening. Fujita 
et al. (2013) presented a scenario that cosmic ray accel¬ 
eration could consume the kinetic energy of infalling gas 
and affect the entropy profile in the cluster outskirts. An¬ 
other explanation is provided by Cavaliere et al. (2011), 
who propose a cluster evolutionary model incorporating 
the effects of merger shocks weakening over time (towards 
low redshift); Walker et al. (2012a) finds these models 
are consistent with the observed entropy flattening, al¬ 
though their model contains several free parameters. All 
the above explanations predict that the entropy profile 
depends on the relaxation state and mass of the cluster. 
Most studies have focused on massive, cool-core clusters. 
It is important to extend these studies to lower mass 
clusters. 

Indeed, an intriguing example of such a system is the 
poor-cluster/fossil-group RXJ1159+5531 (Humphrey et 
al. 2012). We obtained good constraints on its gas prop¬ 
erties all the way out to i?ios, finding that at its virial 
entropy profile agrees with the prediction from gravity- 
only simulations and its baryon fraction within Rios is 
fully consistent with the cosmic value (0.15 Planck Col¬ 
laboration 2013; 0.17 Komatsu et al. 2011). Fossil groups 
are empirically defined as systems with 1) a central dom¬ 
inant galaxy more than two optical magnitudes brighter 
than the second brightest galaxy within half a virial ra¬ 
dius; and 2) an extended thermal X-ray halo with T(x,boi) 
> 10 42 hi q 2 erg s -1 (Jones et al. 2003). Perhaps as a 
fossil-group, RXJ1159+5531 may be sufficiently evolved 
and relaxed (e.g. Harrison et al. 2012) that hydrostatic 
equilibrium is an accurate approximation, and there is 
little clumping. We note that the entropy profile of 
ES03060170, the only other fossil group that has been 
observed with Suzaku out to R, v w, does flatten, but its 
deviation occurs at a much larger radius compared with 
other systems, and the value of the entropy at i ?200 is still 
consistent with the predictions from gravity-only simu¬ 
lations (Su et al. 2013). 

Recent Suzaku observations of cluster outskirts indi¬ 
cate that azimuthal asymmetries are common, even in 
many clusters that appear symmetric at small scales 
(e.g., Eckert et al. 2013; Urban et al. 2014). If 
RXJ1159+5531 is a highly evolved and relaxed system, 


we would, in contrast, expect to see little azimuthal 
variation in these properties as predicted in simulations 
(Vazza et al. 2011). The Suzaku observation offset to the 
North reported in Humphrey et al. (2012) only provides 
~ 27% azimuthal coverage beyond i?soo- To achieve a 
complete azimuthal coverage of this valuable system, we 
acquired deep Suzaku observations in the other three di¬ 
rections. Together, these observations allow the entire 
virial radius to be studied. 

In the cluster outskirts, where the X-ray emission is 
background dominated, the low and stable instrumen¬ 
tal background of Suzaku is crucial for constraining the 
properties of the hot gas. However, the observed flat en¬ 
tropy profile if, due to gas clumping, would result from 
substructures unresolved by the point-spread function 
(PSF) of Suzaku. As an independent study, Eckert et 
al. (2013) combined Planck pressure profile and ROSAT 
density profile and obtained entropy profiles of cool-core 
clusters in line with the baseline entropy profile. More¬ 
over, in this regime, the cosmic X-ray background (CXB) 
dominates the cluster outskirts, in particular for energies 
above about 2keV. An accurate characterization of the 
background is therefore a prerequisite for a reliable mea¬ 
surement of the gas properties. 

The Chandra X-ray Observatory has superb spatial 
resolution (0.5"), which makes it ideal to address these 
issues and to complement the Suzaku data. In addition, 
the large field-of-view of the Chandra ACIS-I combined 
with a low background (compared to ACIS-S) make it 
very suitable to investigate the cluster outskirts and pro¬ 
vide an independent test of the Suzaku results. A grow¬ 
ing number of Suzaku observations have been awarded 
Chandra follow up observations and have been used to 
refine the analysis of Suzaku data. Miller et al. (2012) 
demonstrated that a short snapshot Chandra observation 
is able to reduce the uncertainties in the surface bright¬ 
ness measured by Suzaku by 50% at the outskirts of clus¬ 
ters. We acquired deep ACIS-I observations, mosaicing 
the portion of the i? 5 oo _ ^vir region in the north direc¬ 
tion that was observed and reported in Humphrey et al. 
(2012). These deep exposures allow us to obtain a better 
characterization of the CXB and allow us to reduce the 
uncertainties of the gas properties in the spectroscopic 
analysis directly. 

Adopting a redshift of z = 0.081 from the NASA/IPAC 
Extragalactic Database (NED), we derive a luminosity 
distance of 368 Mpc (so 1' = 90 kpc), assuming a cos¬ 
mology with H 0 = 70 km s -1 Mpc -1 , Ha = 0.7 and 
fl m = 0.3. We studied this system out to Rios = 12' 
(1100 kpc) in the north, south, west, and east directions 
with Suzaku (R 50 o = 580kpc and i? 2 oo = 871 kpc). We 
determined its gas properties at large radii in all direc¬ 
tions, which we present in this paper. We will present 
the hydrostatic constraints of its total mass distribution 
and dark matter properties in paper II and its metallic- 
ity distributions in paper III. We describe the observa¬ 
tions and data reduction in §2 and introduce our mass 
modeling techniques in §3. We report results in §4 and 
construct detailed systematic error budgets in §5. The 
implications of our results are discussed in §6, and our 
main conclusions are summarized in §7. Uncertainties 
reported in this paper are at ler confidence level unless 
stated otherwise. 
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Fig. 1.— Mosaic of the Suzaku XIS1 images of RXJ1159+5531 
in the 0.5—4.0 keV energy range with NXB subtracted. The image 
is corrected for exposure and vignetting. A region of enhanced 
brightness to the east of the center that contains several bright 
point sources when compared to the Chandra image. Point sources 
marked in yellow circles were excluded. The blue circle indicates 

Rios- 

2. OBSERVATIONS AND DATA REDUCTION 

RXJ1159+5531 has been observed with Suzaku from 
the center and out to the R v i r to north, south, west, and 
east directions. A mosaic of Suzaku pointings is shown 
in Figure [T] RXJ1159+5531 has been observed with 
Chandra with one ACIS-S pointing at the center and two 
ACIS-I offset pointings covering the entire field-of-view 
of the Suzaku north pointi ng 3 4 ! A mosaic of Chandra 
pointings is shown in Figure [2] The result of the Suzaku 
analysis of the north pointing combined with the Chan¬ 
dra analysis of the central region has been presented in 
Humphrey et al. (2012). We processed (and also repro¬ 
cessed) all these Suzaku and Chandra observations to 
guarantee the latest calibrations and consistent reduc¬ 
tion process. The observation logs are listed in Table 1. 
Basically, gas properties beyond 0.5 R v j r of the north di¬ 
rection are obtained with a joint Suzaku and Chandra 
analysis; that of the other three directions were obtained 
with Suzaku observations only. 

2.1. Chandra 

We used the CIAO 4.6 and Heasoft 6.15 software, and 
the Chandra calibration database ( Caldb ) 4.6.3, to re¬ 
duce the data. All data were reprocessed from level 1 
events following the standard data reduction thread^ 
Light curves were extracted from a low surface bright¬ 
ness region of the CCDs using the CIAO script lc_clean. 
High background intervals exceeding 2 a above the mean 
quiescent background rates were excised. Effective expo¬ 
sure times are listed in Table 1. 

3 There is another 19 ksec ACIS-I pointing at the center of 
RXJ1159+5531. We did not include it in our analysis in order to 
to simplify the analysis (in particular, the background modeling) 
and avoid additional systematic uncertainty. 

4 http://cxc.harvard.edu/ciao/threads/index.html 
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Fig. 2.— Mosaic of the Chandra images of RXJ1159+5531 in 
0.3-7.0 keV band corrected for exposure variations, containing one 
central ACIS-S pointing and two ACIS-I offset pointings to the 
north direction. The blue circle indicates -Rios- 

Point sources were detected in a 0.3-7.0 keV image 
with wavdetect, supplied with a 1.7 keV exposure map 
to minimize spurious detections at chip boundaries. The 
detection threshold was set to 10 -6 and the scales of 
wavdetect were set to a +2 series from 1 to 8. All de¬ 
tected point sources were inspected by eye, and corre¬ 
sponding elliptical regions containing 99% of the source 
photons were generated. We obtained the 2.0-8.0 keV 
count rates for each resolved point source. An elliptical 
annulus centered on each point source was used for “lo¬ 
cal” background subtraction, which is just outside the 
source-extraction region. We further converted the net 
count rate into flux for each point source assuming a 
power law spectrum with an index of 1.41 (De Luca & 
Molendi 2004). 

We extracted spectra from seven continuous annular 
regions from the central Chandra observation on the 
ACIS-S3 chip, 0-7", 7-15", 15-25", 25-42", 42-78", 78- 
138", 138-208", and 208-300". The width of each annu¬ 
lus was chosen to contain sufficient photons for spectral 
analysis with approximately the same number (> 2000) 
of background-subtracted counts. Spectral response ma¬ 
trices were produced for each annulus with the CIAO tool 
mkwarf and mkacisrmf. Spectral fitting was performed 
with XSPEC 12.7 using the C-statistic. We rebinned the 
spectra to ensure at least 20 photons per bin to aid in 
model selection and computational speed. The energy 
range for spectral fitting was restricted to 0.5-7.0 keV. 
All spectra were fitted simultaneously including compo¬ 
nents for the cluster emission and background. We used a 
single thermal vapec component to model the ICM emis¬ 
sion. We included an additional 7.3 keV bremsstrauhlung 
component to account for the emission from unresolved 
low-mass X-ray binaries (LMXBs) in the central galaxy 
(Irwin et al. 2003). Since the number of LMXBs scales 
with the stellar light, the relative normalization of this 
component between each annulus was fixed to match the 
relative AT-band luminosity in the associated regions. 

To account for the X-ray background, we employed 
a multi-component background model consisting of an 
apec thermal emission model for the Local Bubble 
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TABLE 1 
Observation log 


Name 

Obs ID 

Obs Date 

Exposure 

(ksec) 

R.A. 

Dec. 

Dectectors 

Suzaku N 

804051010 

2009-05-02 

84 

11 59 48.72 

55 36 39.6 

XIS0,1,3 

Suzaku S 

807064010 

2012-05-27 

81 

11 59 51.29 

55 24 44.3 

XIS0,1,3 

807064020 

2012-12-18 

21 

11 59 55.49 

55 26 24.7 

XIS0,1,3 

Suzaku E 

809063010 

2014-05-29 

96 

12 00 46.70 

55 31 16.0 

XIS0,1,3 

Suzaku W 

809064010 

2014-05-31 

94 

11 58 59.83 

55 32 39.1 

XIS0,1,3 

Chandra Cen 

4964 

2004-02-11 

76 

11 59 51.40 

55 32 01.0 

ACIS-S 

Chandra NW 

14026 

2012-08-09 

50 

12 00 41.60 

55 39 55.1 

ACIS-I 

Chandra NE 

14473 

2012-08-12 

37 

11 59 29.10 

55 42 13.7 

ACIS-I 


14027 

2012-08-09 

13 



ACIS-I 


(apecLB, kT=0.08 keV, solar abundances), an additional 
apec thermal emission model (apecywi kT=0.2 keV, so¬ 
lar abundances) for the Milky Way emission in the line 
of sight (Smith et al. 2001), and a power law model 
powcxB (with index F = 1.41) characterizing the unre¬ 
solved CXB (De Luca & Molendi 2004). All these com¬ 
ponents but the Local Bubble were assumed to be ab¬ 
sorbed by foreground (Galactic) cooler gas, with the ab¬ 
sorption characterized by the phabs model for photoelec¬ 
tric absorption. Photoionization cross-sections were from 
Balucinska-Church & McCammon (1992). We adopted 
a Galactic hydrogen column of Nr = 1-02 x 10 20 cm -2 
toward RXJ1159+5531, deduced from the Dickey and 
Lockman (1990) map incorporated in the HEASARC Nr 
tool. The ICM components were allowed to vary inde¬ 
pendently for each radial annulus. The normalization 
of each background component within each annulus was 
linked to scale with the extraction area, but the total 
normalizations were allowed to vary freely. Best-fit X- 
ray background components are listed in Table 2. To 
accommodate the particle background, we included a 
number of Gaussian lines and a broken power-law model 
(see Humphrey et al. 2012 for detail), which were not 
folded through the ancillary response file (ARF). The 
normalization and shape of the instrumental components 
were allowed to vary freely. 

2.2. Suzaku 

Suzaku data reduction and analysis were per¬ 
formed with the Heasoft 6.15 software package using 
CalDB20141001. Data were obtained in both 3x3 and 
5x5 data readout modes; 5x5 mode data were converted 
to 3 x 3 mode and merged with the 3x3 mode data. The 
events were filtered by retaining those with a geomag¬ 
netic cutoff rigidity > 6 GeV/c, and an Earth elevation 
> 10°. The calibration source regions and hot pixels 
were excluded. Light curves were filtered using CL404.6 
script lc_clean. No anomalous event rates deviate more 
than 3tr from the mean was found. The effective exposure 
time of each pointing is listed in Table 1. Bright point 
sources were identified by eye and confirmed by Chan¬ 
dra imaging whenever possible. We give the number of 
each point source as labeled in Figure |T| We excluded a 
circular region of 2.5' radius centered on each source (we 
excluded a elliptical region of similar size for source #1). 

We extracted spectra from four annular regions cen¬ 
tered on RXJ1159+5531 and extending out to R v j r : 0'- 
2', 2! -5', 5'-9', and 9'-13' from the Suzaku observation in 


our previous analysis of the North pointing (Humphrey 
et al. 2012). In this work, we employed the same radii 
for the North and West pointings. However, we were un¬ 
able to obtain reliable constraints on the ICM properties 
using the same radii for the South and East pointings 
due to their farther offset and point source contamina¬ 
tions. For the South pointing, we used spectra extracted 
from O'—2', 2'-4', 4'-7', and 7'-ll' annulus regions to 
model the ICM properties (and we used its 11'—15' an¬ 
nulus spectrum to assist in constraining the background 
model for this pointing). For the East pointing, we find 
it is necessary to use a smaller number of radius sets of 
larger extraction regions: 0'-2.5', 2.5'- 7', and 7'-13' to 
model the ICM. The FTOOL xissimarfgen was used 
to generate an ARF for each region and detector. To 
provide the appropriate photon weighting for each ARF, 
we used a /3-model surface brightness distribution deter¬ 
mined by fitting the surface brightness profile from the 
central Chandra observations. To model the X-ray back¬ 
ground, we followed the standard procedure and gener¬ 
ated ARFs assuming uniform sky emission with a radius 
of 20'. Redistribution matrix files (RMFs) were gen¬ 
erated for each region and detector using the FTOOL 
xisrmfgen. Non-X-ray background (NXB) spectra were 
generated with the FTOOL xisnxbgen. Spectra from 
XIS0, 1, & 3 were simultaneously fitted with XSPEC 
vl2.7.2 (Arnaud 1996). We adopted the solar abundance 
standard of Asplund at al. (2006) in thermal spectral 
models. Energy bands were restricted to 0.5-7.0 keV for 
the back-illuminated CCD (XIS1) and 0.6-7.0 keV for 
the front-illuminated CCDs (XIS0, XIS3), where the re¬ 
sponses are best calibrated (Mitsuda et al. 2007). In 
order to account for spectral mixing between each annu¬ 
lus, we produced “crossing” ARFs between each radial 
annulus using the algorithm described in Humphrey et 
al. (2011). 

We used XSPEC to fit each spectrum with a multi¬ 
component background plus source model (as for the 
Chandra data). Since the Chandra analysis reveals a 
significant cool core, we included a second apec compo¬ 
nent to better model the ICM emission in the central bin. 
We were unable to constrain the ICM metal abundance 
of the outer most bin in each direction separately; we 
linked the metal abundance of the outer most bin to the 
adjacent annulus interior to it for each direction. To ac¬ 
count for the discrepancies in the responses between the 
front-illuminated CCDs and back-illuminated CCDs, the 
normalizations of the model for Bi-chip (XIS1) and FI- 
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TABLE 2 

X-ray Background Components* 



CXB§ 

Local Bubble^ 

Milky Way* 

Chandra-Center 

3 9 +1 ’ u 

°- y -0.4 


8'9-2.7 

Suzaku-North 

3 g+0- 2 

!6-7 1 2 2 5 3 


Suzaku-South 

8 4+0-3 
°- 4 -0.3 

44.8±ll 

10 -5±i:! 

Suzaku-East 

q 3+0.2 
—0.3 

62 . 2+11 

12-9±?i 

Suzaku-West 

7 9+0-3 

' - y -0.4 

58. 4 13-5 

i2-i ±1;| 


*: Results for the normalizations of the different components of the (non-instrumental) X-ray background for each observation. 

§: Normalization of a power-law component with fixed slope (T=1.4) divided by the solid angle, in the unit of photons s —1 cm —2 keV —1 str —1 
at IkeV. Note that the result for the Suzaku North direction is the remaining unresolved flux after accounting for the sources resolved by 
Chandra. 

An unabsorbed apec thermal component (kT=0.08 keV, solar abundances) with normalization expressed as an emission measure 
integrated over the line of sight, 4 7r [.D A (i-)- z )]2 I n enndV/dfl in the unit of 10 —14 cm -5 str -1 , where dfl is the solid angle. 

An absorbed apec thermal component (kT=0.2 keV, solar abundances) with normalization expressed as an emission measure integrated 
over the line of sight, 4 7r [.p A (i-|-,z )]2 / n e nndV/dQ in the unit of 10 —14 cm -5 str — 1 , where dfl is the solid angle. 



East 7.5-13’ 


West 9-13’ 




Fig. 3.— Suzaku XIS1 spectra for RXJ1159+5531 from the outermost bin of the north, south, west, and east directions. Instrumental 
background has been subtracted. Black: observed spectra. Red: best-fit model. Green: the X-ray background. Magenta, Blue, Yellow, 
Orange: hot gas emission from 1st, 2nd, 3rd, and 4th radial annulus, respectively, [see the electronic edition of the journal for a color 
version of this figure.] 
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chips (XISO, XIS3) were allowed to vary independently; 
the discrepancy of their normalizations is ~10%. In Fig¬ 
ure [3j we show the XIS1 spectra of the outmost annulus 
in each direction and their individual components. 


2.3. Refine Suzaku analysis with point sources resolved 
by Chandra 

The unresolved CXB component is the dominant and 
most uncertain component of the X-ray background in 
the Suzaku spectra of cluster outskirts for energies above 
1 keV. In order to constrain this component, we followed 
the approach described in Walker et al. (2013). We ob¬ 
tained the coordinates and 2.0-10.0 keV count rates of 
the point sources resolved by the three Chandra point¬ 
ings (Figure [2j . We convert their count rates to fluxes 
in 2.0-10.0 kew assuming a absorbed powerlaw spectrum 
with an index of 1.41. We detected 148 point sources in 
the FOV of Chandra and 85 of them are within the FOV 
of Suzaku. The faintest resolved point source has a flux 
of 4.01 xlO -17 ergs -1 cm -2 . 

The surface brightness profile of point sources in 
Suzaku FOV is shown in Figure [4] Using the Heasoft 
tool xissim, we simulated a Suzaku observation of point 
sources resolved with Chandra in the field-of-view of the 
Suzaku observation. The simulated exposure time is set 
to a very large number to guarantee good statistics. We 
extracted spectra from each annular section used in the 
Suzaku analysis from this simulated observation and fit 
the spectra to a absorbed powerlaw with an index of 
1.41, we obtained the total flux of these simulated point 
sources within each annulus. Their corresponding sur¬ 
face brightness profiles are shown in Figure [4]4eft. 

In our Suzaku spectral fitting of the north direction, 
we include a powerlaw component with an index of 1.41 
to account for the emission of the Chandra resolved point 
sources and with their normalizations fixed at their as¬ 
sociated fluxes. We also included another powerlaw 
component with an index of 1.41 to model the remain¬ 
ing unresolved point sources; its normalization is set to 
vary freely but all regions were linked together since the 
fluctuations of the remaining unresolved point sources 
should be much smaller. We obtained a best-fit of 
7.14x 10 -12 ergs -1 cm' 2 deg -2 for the remaining unre¬ 
solved point sources in 2.0-10.0 keV. 

Following Moretti et al. (2003), the expected level of 
the remaining unresolved point sources can be estimated 
through the integral, 

/ c /]y\ 

Fcxb = 2.18 ±0.13 xlO" 11 - / 1—jxSdS, (1) 

in units of ergs -1 cm -2 deg -2 . We take an analytical 
form of the source flux distribution in the 2.0-10.0 keV 
band given by Moretti et al. (2003), 

N{> S) = 


Nsm _ < 2 * IQ- 18 )— _ (2) 

where a hH = l-57±g;Jg, a 2M = 0.44^;^, S 0 , H = 
(4.5+ 3 ; 7 ) x 10 -15 ergs -1 cm -2 , and Nh = 5300^400. 
Using the best-fit values of these parameters, we esti¬ 
mate the level of remaining unresolved point sources to 


be 6.69 ± 3.76 x 10 -12 ergs -1 cm -2 deg -2 . The best- 
fit surface brightness of unresolved point sources are in 
good agreement with the expected emission of unresolved 
point sources as compared in Figure [4]4eft. 

For most annuli, about 50% of the point sources were 
resolved (see Figure pjdeft). The total flux of the resolved 
point sources plus tie flux of the remaining unresolved 
point sources are comparable to the flux for the originally 
determined unresolved point sources with Suzaku data 
only. In Figure [J]-right we compare the Suzaku measure¬ 
ments of the entropy profile of RXJ1159±5531 (see below 
in §4.3) in the north direction, obtained with and with¬ 
out the refinement of Chandra. The values are consistent 
while the statistical uncertainty of the virial entropy has 
been reduced by 18% with the incorporation of Chandra 
data. 


3. HYDROSTATIC MODELS 

The forward fitting approach allows us to fit the pro¬ 
jected density and temperature profiles directly and also 
allows the cluster emission outside the largest spectral 
extraction annulus to be treated self-consistently. These 
properties are advantageous for background-dominated 
X-ray emission in the outskirts of a cluster where spec¬ 
tral deprojection (e.g., using the standard onion-peeling 
approach) further degrades the data quality and tradi¬ 
tionally assumes there is no cluster emission outside the 
largest annulus used. We focus on solutions of the hy¬ 
drostatic equilibrium equation in terms of the entropy 
because the additional constraint of convective stability 
is easily enforced by requiring the entropy to increase 
with radius. In addition, unlike the temperature and 
density, the much smoother entropy profile is easily pa¬ 
rameterized with a simple broken power law model. 

Rewriting the hydrostatic equation in terms of the en¬ 
tropy and a variable depending on the pressure yields 

^ = - - GM( < r) K~ 3 ' 5 . £ = P 2/5 (3) 

dr 5 r z 

where P is the gas pressure, and M(< r) is the total 
enclosed mass within radius r, which includes the contri¬ 
butions from stars, gas, and dark matter. Equation (4) 
can be solved directly given models for K(r) and M(r) 
provided the gas mass can be neglected, since the pres¬ 
sure (and therefore £) depends on the gas density. For 
a self-consistent solution Equation (4) is rearranged and 
differentiated with respect to r, taking advantage of the 
mass continuity equation to give, 


— f r 2 /\ 3/5 — 


dr 


8tt r 2 G 


dr 
8 tt r 2 G 


k - 3 / 5 ^ 3/2 


(.Ps 


Pdm )• 


( 4 ) 


To parameterize the dark matter component we used 
the NFW profile (Navarro et al. 1997), 


PdmM = 


Po 


i (! + *)' 


( 5 ) 


where the free parameters are the scale radius ( R s ) and 
a characteristic density (po)- The stellar component has 
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Fig. 4. — left: black: surface brightness profile of point sources resolved by Chandra, red: surface brightness profile of simulated point 
sources for Suzaku. Green: best-fit surface brightness profile of the CXB component determined with Suzaku. Magenta: best-fit surface 
brightness profile of the remaining CXB component determined by spectral fitting with Suzaku after excluding point sources resolved by 
Chandra. Blue: the expected surface brightness of the remaining CXB component (dashed blue: lcr error), right: Entropy profile of 
the north direction. Blue: results measured with Suzaku observation only. Green: results measured with Suzaku and the point source 
refinement of Chandra. 


no free parameters and is described in §4.4. For K(r) we 
adopt a simple broken power-law model plus a constant, 
which takes the form of K = ko + /(r), where 


ki r^ 1 

r < n 

k /2 r@ 2 

7"1 < r < 7~2 

k 3 r& 3 

r > r 2 


There are two boundary conditions that need to be spec¬ 
ified for £(r) and d£(r)/dr. Since the gas mass is negligi¬ 
ble at small radius, Equation (4) can be solved directly 
at some small radius to give the boundary condition for 
dt;(r)/dr. The boundary condition for £(r) amounts to 
specifying the pressure or entropy at some radius which 
we treat as a free parameter (see Humphrey et al. 2008 
and Buote & Humphrey 2012). For a given set of param¬ 
eters for K(r), M(r), and the £(r) boundary condition, 
we solve Equation (7) for £(r) from which we compute the 
density (emission measure) and emission-weighted tem¬ 
perature profiles in projection that are then compared to 
the measurements via the y 2 statistic. We explored the 
parameter space using a Bayesian Monte Carlo Method 
— version 2.7 of the MultiNest code (Feroz & Hobson 
2008; Feroz et al. 2009). Flat priors were assumed for all 
free parameters. 

In this work, we applied a pair of parameterized dark 
matter density and entropy profiles that can uniquely 
determine the pressure profile by solving the hydrostatic 
equilibrium equation (there are also free parameters as¬ 
sociated with the boundary conditions). We can ob¬ 
tain three dimensional temperature and density profiles 
through entropy and pressure profiles. We project the 
three dimensional temperature and density profiles on to 
the sky assuming spherical symmetry of the ICM and 
compare them to the observed projected temperature 
and density profiles. We repeat this process through 
monte-carlo simulations until we can reproduce the ob¬ 
served results. The ability to fit the projected data di¬ 


rectly is very valuable for low quality, noisy data such 
as cluster outskirts whose integrity would be greatly de¬ 
graded by deprojection noise. In this work, we applied 
this technique to data in each direction separately to ob¬ 
tain their gas properties. We display the projected tem¬ 
perature and density profiles along with their best-fitting 
models in Figure [5] (top and bottom panels). We ob¬ 
tained a best-fit y^/d.o.f of 13.5/11, 12.6/11, 10.8/9, and 
12.8/11 for the North, South, East, and West directions 
respectively. As shown in §4, very similar results for the 
entropy and other quantities are obtained for each direc¬ 
tion. Thus, for convenience and clarity of presentation 
we adopt fiducial values of R 500 , i? 200 j -Rios, and the to¬ 
tal mass within i?ios (M to t=9.7x 10 13 M 0 ) obtained from 
the North pointing as the parameters of RXJ1159+5531. 
We note that in our joint fit of the Chandra and Suzaku 
data the high S/N Chandra data at small radii are crucial 
for constraining the model components near the center 
(e.g., M/L of the stellar component, NFW scale radius). 
However, the Suzaku data provide the crucial constraints 
on the model properties obtained at large radius; the 
constraints on the mass profile derived using only the 
Chandra data (e.g., Gastaldello et al. 2007) were not as 
good as our joint fit. 

4. RESULTS 

4.1. Temperature and density profiles 

The best-fit two-dimensional (2D) gas temperature and 
density profiles of RXJ1159+5531 in all directions are 
shown in Figure [5] (top and middle panels). For each 
direction, Suzaku results are obtained with ~ 90° sec¬ 
tional ~3.0' radial bins beyond 0.3 R vlr . We did not 
plot the results of the innermost bin of Suzaku since our 
two-temperature fit is only a crude representation of the 
spatially resolved temperature profile measured by Chan¬ 
dra. We also plotted the results obtained with Chandra 
in several ~ 1', ~ 360° annuli for the central regions. 
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Fig. 5. — Top: Projected density profiles of central regions measured with Chandra (black diamonds) and of each direction measured with 
Suzaku (green diamonds). Black and green solid lines indicate the values of best-fits. Bottom: Projected temperature profiles of central 
regions measured with Chandra (red) and of each direction measured with Suzaku (blue), red and blue solid lines indicate the values of 
best-fits. 


The combined profiles span from the center out to R v i r . 
The width of our extracted annuli vary from 4 kpc in 
the center to 350 kpc out to i? v ir- In the overlapping re¬ 
gions, the results from Chandra and Suzaku show good 
agreement. 

We plot the three-dimensional (3D) density and tem¬ 
perature profiles for each direction in Figure [6] and Fig¬ 
ure [d respectively (marked in solid lines) obtained with 
the forward-fitting approach described in §3. We com¬ 
pare the density and temperature profiles of all directions 
together in Figure [8 (a) and (b). Each direction shows 
remarkably similar aensity profiles. The virial temper¬ 
atures of all directions are consistent within lcr uncer¬ 
tainties. The 3D temperature profiles decline by more 
than a factor of 2 from a peak temperature of ^3kcV 
at 0.3 R vir to ~1.0keV at R v i r . We compared the 3D 
temperature profiles to the heuristic formula calibrated 
by Pratt et al. (2007) for clusters over a radial range of 
0.125 < R/R 200 < 0.5: 


ployed a power-law component with two breaks and a 
constant component. The normalizations of the power- 
law and constant components, the radii of the breaks 
and their slopes were all free parameters. We plot the 
best-fit model of the 3D entropy profile for each direc¬ 
tion in Figure [9] Shaded regions indicate lcr confidence 
regions. We marked in magenta the baseline predictions 
from gravitational structure formation (Voit et al. 2005), 
given by: 

Agra(i?.) = 1.32 K 2O0 (R/R 2 oo) 11 , (8) 

where the normalization, K 2 oo^ is given by 


-ft 200 — 362 


GM 2 ooiJ,m p 

2i?200 



H(z) 

Hn 


1-4/3 


03 


- 4/3 


keV cm 


-2 


T/T x = 1.19 - 0.74 R/R 200 , (7) 

where Tx is the peak temperature. We observe a similar 
decline of the temperature profiles out to large radii in 
this system as observed in other galaxy clusters (Aka- 
rnatsu et al. 2011). 

4.2. Entropy and pressure profiles 

We fitted the broken power-law model (equation 6) for 
the entropy to the data and enforced a monotonically 
rising profile (for convective stability) by requiring the 
slope parameters to be positive or zero. This model em- 


The central entropy profile is more elevated and extended 
than this baseline A' gra . The entropy profiles of the north 
and east directions beyond 0.5 R v i r are very well behaved, 
rising linearly all the way out to the R v ; r , following the 
K oc r 11 expectation. The entropy profiles of the south 
and west directions have a flatter slope but stay consis¬ 
tent wit h th e theoretical expectation near i? v i r - 
Figure [To] shows the three dimensional pressure profile 
of each direction. We compared the observed pressure 
profile of each direction to a semi-analytic universal pres¬ 
sure profile derived by Arnaud et al. (2010) from com¬ 
parison of their numerical simulations to XMM-Newton 
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Fig. 6.— Solid lines: 3D density profiles derived with our primary method (entropy-based forward fitting) for each direction. Shaded 
regions indicate lcr uncertainties. Triangles: 3D density profiles obtained via spectral deprojection pro jet (Brown measured with Chandra ; 
Black measured with Suzaku). 


observations of clusters within P 50 o. This pressure pro¬ 
file is characterized as 


P(r) = Pi 


500 


M; 


500 


-| a-p+a' (x) 


1.3 X 10 14 /l7Q M 0 
Po 


(c 50 oa;) 7 [l + (c50oa;)“] (/3 “ 7)/a ’ 


where 


(9) 


< (l ) =0 - 10 - 1 °’’ +0 - 1 °) [ 1 +('°/o ) 5)»r 


and P 500 and M 500 are respectively the pressure and total 
mass at P 50 o. Arnaud et al. (2010) adopted parameter 
values of 

[Po, c 500 , 7, P\ = 

[8.403/i7 0 3/2 , 1.177, 0.3081, 1.0510, 5.4905]. 

The pressure profiles of the all direction of 
RXJ1159+5531 from 100 kpc out to about half P 500 
are in good agreement with this universal profile. How¬ 
ever, its pressure profiles exceed this universal profile at 
P 500 by 60% to 200 %. Note that this universal pro¬ 
file is an average result of a large number of different 
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Fig. 7. — Solid lines: 3D temperature profiles derived with our primary method (entropy-based forward fitting) for each direction. Shaded 
regions indicate 1 a uncertainties. Magenta: the heuristic formula calibrated by Pratt et al. (2007) for clusters. Triangles: 3D temperature 
profiles obtained via spectral deprojection pro jet (Brown measured with Chandra; Black measured with Suzaku). 


clusters which have a 300% internal discrepancy at -R 500 
itself. We do not expect the observed pressure profiles 
of RXJ1159+5531 to be in perfect agreement with this 
average profile. At least, this result suggests that the 
pressure profiles of RXJ1159+5531 are not significantly 
smaller than other systems and its outskirts are unlikely 
to be dominated by non-thermal pressure support. We 
compared the pressure profile of all directions in Figure [8] 
(d). They have remarkably similar behaviors out to the 
large radii. 

4.3. Gas and baryon fractions 


To compute the gas and baryon mass, we included the 
stellar mass of the central galaxy, the gas mass, and the 
two additional baryon reservoirs: intracluster light and 
other member galaxies. Vikhlinin et al. (1999) found 
that ~25% of the F-band stellar light is associated with 
other member galaxies. We assumed a similar ratio in 
the A'-band and adopted a M*/A K ratio of 1 for these 
galaxies. Furthermore, we assume that the intracluster 
light within i? v j r contains twice as much as the stellar 
mass of the central galaxy (Purcell et al. 2007). 

A 200 (871 kpc) is near the effective center of the out¬ 
ermost radial bin in all directions. We confirm that we 
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Fig. 8. — Gas properties of all directions obtained with our primary method (entropy-based forward fitting). Shaded regions indicate lcr 
uncertainties, (a). 3D density profile, (b). 3D temperature profile, (c). 3D entropy profile. Magenta: entropy profile from gravity-only 
cosmology simulations, K oc r 1 ' 1 ; normalization is derived from Voit et al. (2005). (d). 3D pressure profile. Magenta: universal pressure 
profile derived by Arnaud et al. (2010). 


are able to measure the gas properties at i ?200 with these 
data as listed in Table 3. The baryon fractions of each di¬ 
rection are consistent within uncertainties. i?ios is within 
the outermost bins of north, east, and west direction, 
and just outside that of the south direction. Based on 
our entropy-based forward-fitting method, we can put 
constrain on its gas properties around i?ios- The radial 
profiles of enclosed gas mass fraction and baryon fraction 
of each direction are shown in Figure [TTI We obtained 
total masses within i?ios of 9.72±1.44 (N), 9.91±1.74 
(S), 9.66T2.10 (E), and 10.9T2.0 (W)xlO 13 M 0 for each 


direction as described in §3. The enclosed baryon frac¬ 
tions within R W8 are O.lGllg^g (N), 0.155 +q;q|| (S), 

0.162+q° 44 (E), and 0.158jjgQ37 (W) for each direction. 
All of them are very consistent with the cosmic baryon 
fraction (0.15 Plank Collaboration 2013; 0.17 Komatsu 
et al. 2011). 

5. SYSTEMATIC ERROR BUDGET 

We have constructed a detailed error budget consid¬ 
ering a variety of possible systematic effects. In Ta¬ 
ble 3, we list error budgets for the densities, tempera¬ 
tures, entropies and pressures measured at R 200 , as well 
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Fig. 9.— Solid lines: 3D entropy profiles derived with our primary method (entropy-based forward fitting) for each direction. Shaded 
regions indicate la uncertainties. Magenta line: entropy profile from simulations, S oc r 11 ; normalization is derived from Voit et al. (2005). 
Triangles: 3D density profiles obtained via spectral deprojection pro jet (Brown measured with Chandra; Black measured with Suzaku). 


as enclosed baryon and gas fractions within R 200 , for 
the north, south, west, and east directions. We consider 
any systematic error to be significant once it equals or 
exceeds the statistical error of the same variable. Our 
discussion pays special attention to the impact of each 
effect on the measurement of the entropy at i? 200 - 

5.1. Background 

In the outermost annulus, the background dominates 
the X-ray emission so that the ICM contributes only 8 
(1)%, 11 (1)%, 35 (0.2)%, and 6 (0.4)% of the total emis¬ 
sion (including the NXB) in 0.5-2.0 keV (2.0-8.0 keV) 
for north, south, east, and west directions, respectively. 


Consequently, our results depend on the determination 
of the various background components. To examine our 
sensitivity to the background, we first increased and de¬ 
creased the level of particle background by 5% since the 
variation of the XIS NXB is expected to be <3% (Tawa 
et al. 2008). We list the impact of these variations in Ta¬ 
ble 3 (as ANXB). The systematic uncertainty associated 
with variations in the particle background are smaller 
than the statistical uncertainty on the parameters of all 
directions. 

For Suzaku, the most important and also the most un¬ 
certain component of the X-ray background is the CXB. 
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Fig. 10.— Solid lines: 3D pressure profiles derived with our primary method (entropy-based forward fitting) for each direction. Shaded 
regions indicate lcr uncertainties. Magenta: universal pressure profile derived by Arnaud et al. (2010). Triangles: 3D density profiles 
obtained via spectral deprojection pro jet (Brown measured with Chandra; Black measured with Suzaku). 


To assess the systematic uncertainties associated with 
the CXB for each direction, we alternatively fixed the 
normalization of the CXB power law T = 1.41 component 
in the model at the expected value according to Equa¬ 
tion (1) in §2.3 for each direction (see Table 3 ACXB). 
These changes do not have significant impact on our 
results of entropy. Despite that the slope of the CXB 
power law component is well determined (e.g., De Luca 
et al. (2004): T = 1.41 ± 0.06 -XMM-Newton- Moretti 
et al. (2009): T = 1.47 ± 0.07 -Swift; Tozzi et al. (2001): 
T = 1.44± 0.03-Chandra), we examined this uncertainty 
by fixing its slope at T = 1.5 and T = 1.3, respectively. 


The impact of using T = 1.3 appears to be a small effect 
on our results, while using T = 1.5 would increase the en¬ 
tropy at i ?200 by ~ 30% for the south and east pointings 
(see ACXB-T in Table 3). The unresolved CXB com¬ 
ponent may have harder spectrum. Using Chandra and 
Swift, Moretti et al. (2012) determined that the slope 
of the unresolved CXB power law component can be as 
small as 0.1. We tested this effect by applying power law 
T = 0.1 to model the remaining CXB component of the 
north pointing; its results are listed in Table 3 (as ACXB- 
T = 0.1). This variation causes a systematic error on the 
entropy at R 200 slightly larger than its statistical uncer- 
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Fig. 11. — Solid triangles: enclosed gas mass fraction for each direction. Open triangles: enclosed baryon fraction for each direction. 
Horizontal magenta line: the cosmic value determined by Planck. Horizontal black line: the cosmic value determined by WMAP. 


tainty. However, we consider F = 0.1 for the unresolved 
CXB to be an extreme case and its justification requires 
more work in the future. 

5.2. SWCX 

In order to probe the effect of the activity of the So¬ 
lar Wind Charge eXchange (SWCX), we compared the 
light curve of Suzaku observations in the soft band (0.5- 
2.0 keV) to the data taken from the Advanced Com¬ 
position Explorer (ACE) during these observations as 
shown in Figure |12| For the south direction, we man¬ 
aged to obtain the light curve of proton fluxes taken 
from the Solar Wind Electron Proton Alpha Monitor 
(SWEPAM) onboard ACE. Its variation ranges within 


1.5-3.6x 10 8 cm -2 s -1 . According to Snowden et al. 
(2004), a proton flux of < 4x 10 8 cm -2 s _1 is considered 
to be a quiescent level. Unfortunately, the SWEPAM 
data during the Suzaku observations of the other three 
pointings were indicated to be bad or missing. Instead, 
we inspected the 0 +7 /0 +6 ratio taken from the Solar 
Wind Ion Composition Spectrometer (SWICS) onboard 
ACE. Their light curves lack obvious variations and their 
levels (< 0.3) are consistent with the expectations of qui¬ 
escent emission interval (Snowden et al. 2004). Over all, 
the activity of SWCX was inferred to be low during these 
Suzaku observations. Still, we added a few gaussian lines 
to model the SWCX component as part of the X-ray 
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background component, which contains 6 gaussian lines 
at energies of 0.46 keV, 0.56 keV, 0.65 keV, 0.81 keV, 0.91 
keV, and 1.34 keV (Snowden et al. 2004). We compared 
our best-fit when including and excluding this compo¬ 
nent in the model in Table 3. The differences are typ¬ 
ically much smaller than the statistical uncertainty on 
those parameters. 

5.3. PSF 

The large, energy dependent PSF of Suzaku causes 
some photons from the central part of clusters to be scat¬ 
tered by the optics out to radii of ~ 10'. In our analysis, 
we accounted for this effect by including spectral mixing 
between each annulus in the model. To explore how sen¬ 
sitive our results are to the mixing level, we experimented 
with adjusting the amount of light that is scattered into 
each annulus by ±5%. This did not appreciably affect 
the results of the north and west directions, but it did 
have significant impact on the east and south directions 
(change the entropy at -R 200 by ±35% and ±38% respec¬ 
tively), as shown in Table 3. 

5.4. Entropy model 

The entropy profiles we calculated assume a broken 
power-law model consisting of two breaks and a con¬ 
stant. This model has more freedom than the one used in 
Humphrey et al. (2012) which consists of only one break. 
In order to test how sensitive our result is to the chosen 
entropy model, we added one more break in the model 
of entropy profile. We listed the corresponding results 
in Table 3 (Amodel). This effect causes little impact on 
our results of the east, south, and west directions but it 
gives a higher entropy value for the north direction. 

5.5. Solar abundance table 

Plasma emission is a strong function of metallicity for 
low temperature systems (particularly for kT < 1 keV). 
Adopting an accurate solar abundance table is crucial for 
determining both the thermal properties and the metal 
abundance of hot gas (e.g., Su & Irwin 2013). The dis¬ 
crepancy in measurements, especially in thermal prop¬ 
erties, caused by using different solar abundance tables 
have been largely underappreciated. Comparing the sys¬ 
tematic uncertainty caused by using different solar abun¬ 
dance tables allow us to compare our results to previous 
studies. In our analysis, we adopt Asplund et al. (2006) 
as the abundance table for the plasma emission model. 
We performed our fit again using one of the most out¬ 
dated but widely used solar abundance table Anders & 
Grevesse (1989). We compared these two sets of results 
in Table 3 (Asolar). The discrepancy in the entropy val¬ 
ues is quite large for all pointings. RXJ1159±5531 is 
a low mass cluster for which the contribution of metal 
line emissions to the X-ray emission is more significant 
than more massive clusters. Thus its gas properties are 
more sensitive to the choice of abundance table. Our test 
shows that it is important to use updated solar abun¬ 
dance tables and caution needs to be taken when com¬ 
paring recent studies with previous studies, at least for 
the outermost virial region of low mass clusters where 
kT ~ IkeV. We did the same exercise by using another 
commonly used solar abundance table Lodders (2003); 
the results are also listed in Table 3. Its deviation from 


our best-fit is similar compared with using Anders & 
Grevesse (1989). 

The Fe abundance profile of RXJ1159±5531 of each 
direction are shown in Figure [13| The central dip of the 
metallicity profile observed by TJhandra is due to the “Fe- 
bias” (Buote 2000) in the central region, as has also been 
noted and discussed in Humphrey et al. (2012). We ob¬ 
tained a similar best-fit hot gas metallicity of ~ 0.2 Z 0 
in the spectral analysis for the outermost two annuli in 
each direction (the ICM metallicity of the outermost an¬ 
nulus was tied to their adjacent inner annulus). We ex¬ 
amined the impact of the metallicity determination on 
our results by fixing the hot gas metallicity of the outer¬ 
most two annuli at 0.1 Z 0 and 0.3 Z 0 respectively. Using 
0.3 Z 0 typically has little impact on our results while us¬ 
ing 0.1 Z Q would reduce the entropy of the north and 
south direction by 30%, as shown in Table 3 (Aabun). 

5.6. nH and distance 

To test the sensitivity of our result to the absorbed 
Galactic Ah value. We varied nH from the Dickey and 
Lockman (1990) map by 20%. This did not appreciably 
affect our conclusions except for the east direction as 
shown in Table 3. We also examined the error associated 
with distance uncertainty. We varied the value of redshift 
in the model by 5%. The impact of this latter factor is 
always very small with respect to the statistical errors. 

5.7. FI-BI 

Although there have not been detailed calibration 
works on the level of consistency between different XIS 
chips that we are aware of, we allow the normalizations 
of the XIS1 spectra and that of XIS0 and XIS3 spectra 
to vary independently in the fit in order to account for 
the discrepancies between the response of Fl-chips and 
Bi-chips (we assume the same normalization for the XIS0 
and XIS3). We found a ~10% difference in the normal¬ 
izations determined by Suzaku Fl-chips and Bi-chips for 
RXJ1159±5531. In order to test the significance of this 
discrepancy, we also performed our fit by linking their 
normalizations. The difference it made is smaller than 
the statistical uncertainty in the value of entropy. Details 
are listed in Table 3 (AFI-BI). 

5.8. Techniques 

In order to assess the sensitivity of our results to the 
forward-fitting analysis of the projected temperature and 
density profiles, we compared our results to those ob¬ 
tained using the well-known spectral deprojection model 
pro jet in Xspec. The 3D density and temperature pro¬ 
files obtained with both methods are compared in Fig- 
ures[6]and [7J Using the 3D temperature and density pro¬ 
files derivedwith the pro j ct method, we also calculated 
their corresponding entropy and pressure profiles which 
we also show in Figure [9] and Figure [To] This procedure 
(similar to the classic onion peel) assumes constant ICM 
emission in each three-dimensional shell and no emission 
outside the bounding shell/annulus. 

Both methods have their pros and cons (for a more 
detailed comparison see Buote & Humphrey 2012). A 
key advantage of spectral deprojection with pro j ct is 
that no parameterized model for the radial ICM prop¬ 
erties (e.g., entropy) is required nor is hydrostatic equi¬ 
librium assumed. On the other hand, the forward-fitting 
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Fig. 12.— Crosses: Suzaku XIS1 light curve in 0.5-2.0 keV energy band in each direction, in the unit of counts/s. top-left.: triangles 
represent the light curve of 0 +7 /0 +() ratio taken from ACE-SWICS (1.1) during Suzaku observation of the north direction, top-right: 
circles represent the light curve of proton flux taken from ACE-SWEPAM during Suzaku observation of the south direction, in the unit 
of 10 s cm -2 s _1 . bottom-left: triangles represent the light curve of the 0 +7 /0+ 6 ratio taken from ACE-SWICS (2.0) during Suzaku 
observation of the east direction, bottom-right: same as bottom-left but for the west pointing. ACE data corrected for travel time to Earth, 
[see the electronic edition of the journal for a color version of this figure.] 


method allows us to self-consistently account for ICM 
emission beyond R V1I (we evaluate the hydrostatic model, 
and hence the entropy profile, out to 2f? V i r ) while the 
method using proj ct may over subtract the background 
component by assuming there is no emission outside of 
the outmost extracted region. Moreover, by projecting 
our 3D model onto the sky, the forward-fitting method al¬ 
lows us to correctly account for the radial variation of the 
spectral properties within each bin. Using pro j ct or the 
standard onion-peeling method has to assume the spec¬ 
tral properties are uniform in the entire bin, which can be 


a substantial error for very large apertures such as used 
for Suzaku analysis. We find that the three-dimensional 
radial profiles of the ICM temperature, density, entropy, 
and pressure in each direction obtained by both meth¬ 
ods are consistent with each other as demonstrated in 
Figures [6] [lO] 

6. DISCUSSION 

An elliptical cluster in hydrostatic equilibrium, where 
the ICM is everywhere single-phase, should exhibit no 
substructure in its X-ray emission. While smooth ra- 
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TABLE 3 

Gas properties at R200 with their systematic budgets 




Temperature 

(keV) 

Density 
(10 -5 cm -3 ) 

Entropy 
(keV cm 2 ) 

Pressure 
(10 -5 keV cm -3 ) 

/gas 

/b 


best-fit 

1.21 ± 0.17 

5.98±1.20 

800+145 

7.29+1.90 

fi 12 c i +U ' U1S 
u. izo — 0-01 7 

n 1 QQ+u.uis 

u. ioo_o 018 


ANXB 

±0.06 

±0.60 

±94 

±0.46 

± 0.002 

±0.003 


ACXB 

-0.05 

+0.20 

-55 

+0.12 

+ 0.010 

+0.011 


ACXB-r 

±0.05, -0.10 

-1.18,+0.62 

+135,-115 

-0.94,+0.10 

-0.017,+0.003 

-0.019,+0.003 


ACXB-r (Extreme) 

±0.13 

+4.02 

-180 

+3.01 

+0.023 

+0.024 

North 

ASWCX 

±0.08 

-0.02 

+48 

+0.39 

+0.003 

+0.002 

APSF 

±0.05 

±1.03 

±63 

±1.55 

±0.008 

±0.006 


Amodel 

-0.015 

-0.79 

+62 

-1.09 

-0.001 

-0.003 


Asolar 

-0.012,+0.11 

-0.26,-1.12 

+9.6,+194 

-0.44, -0.87 

-0.005,-0.009 

-0.004, -0.009 


Aabun 

-0.40, -0.01 

-0.20,-0.88 

-256,+75 

-2.60, -1.17 

-0.002,-0.016 

-0.001, -0.017 


AnH 

±0.09 

±0.26 

±86 

±0.30 

±0.04 

±0.03 


ADistance 

±0.01 

±0.15 

±1 

±0.29 

±0.001 

±0 


AFI-BI 

-0.02 

-0.44 

+22 

-0.67 

+0.002 

+0.001 


best-fit 

0.98 ±0.11 

6.77±0.78 

591+90 

6.59+1.1 

U.1ZO_ 0 022 

n 1 00+U.U2U 
U.lOO_ 0 022 


ANXB 

±0.09 

±0.87 

±2.37 

±1.44 

± 0.002 

±0.003 


ACXB 

-0.21 

-1.5 

-43 

-2.62 

-0.006 

-0.008 


ACXB-r 

+0.28,+0.12 

-1.07,+0.57 

+244,+41 

+0.40,+1.48 

-0.017,+0.008 

-0.017,+0.008 


ASWCX 

+0.09 

+0.19 

+67 

+0.34 

+0.005 

+0.004 

South 

APSF 

±0.12 

±1.89 

±227 

±1.2 

±0.015 

±0.015 


Amodel 

-0.01 

+0.66 

-49 

+0.64 

+0.001 

+0.001 


Asolar 

+0.27, -0.01 

-1.38,+0.79 

+287, -50 

+0.02,+0.61 

-0.007,-0.009 

-0.007, -0.009 


Aabun 

-0.27,+0.06 

+0.97,+0.13 

-202,+26 

-1.2,+0.45 

-0.002,-0.017 

-0.002, -0.018 


AnH 

±0.065 

±1.08 

±114 

±0.78 

±0.003 

±0.004 


ADistance 

±0.125 

±1.12 

±159 

±0.49 

±0.003 

±0.003 


AFI-BI 

+0.005 

-0.02 

+2 

-0.07 

+0.002 

+0.002 


best-fit 

1.13 ±0.19 

7.45±1.44 

643+129 

8.43+3.00 

0 126 +UU2y 

U. lZO_Q Q27 

0 i36 _l " u u ^ y 

U. .LOU — 0.027 


ANXB 

±0.11 

±1.39 

±127 

±0.58 

±0.018 

±0.018 


ACXB 

+0.06 

+0.83 

-15 

+1.4 

+0.007 

+0.004 


ACXB-r 

+0.55, -0.01 

+1.51.+2.41 

+187,-110 

+6.35,+2.72 

+0.005,+0.019 

+0.006,+0.015 


ASWCX 

+0.07 

-0.42 

+ 13 

+ 1.00 

0 

-0.002 

East 

APSF 

±0.50 

±0.75 

±223 

±4.94 

±0.012 

±0.013 


Amodel 

0 

-0.64 

-37 

+0.71 

-0.004 

-0.002 


Asolar 

-0.15,+0.08 

+2.14,-1.28 

-174,+134 

+0.97, -0.96 

+0.020,-0.011 

+0.021, -0.012 


Aabun 

-0.06, -0.11 

-0.83,-1.99 

+12,+68 

-1.34, -2.86 

-0.004,-0.023 

-0.005, -0.023 


AnH 

±0.44 

±1.0 

±175 

±4.83 

±0.007 

±0.004 


ADistance 

±0.04 

±0.5 

±52 

±0.24 

±0.005 

±0.004 


AFI-BI 

+0.23 

+1.43 

+42 

+3.65 

+0.020 

+0.020 


best-fit 

1.11 ± 0.12 

7.41+0.74 

634+88 

8.20+1.30 

0 1 9'7+U.U^I 
u - iZ ' -0.021 

0 i34+ u u ^ 

u.io^t_o 022 


ANXB 

±0.005 

±1.51 

±78 

±1.6 

±0.006 

±0.007 


ACXB 

-0.11 

-0.64 

-30 

-1.43 

-0.008 

-0.008 


ACXB-r 

-0.16,+0.03 

-0.65,+1.48 

-56,-51 

-1.81,+1.80 

-0.004,+0.011 

-0.003,+0.012 


ASWCX 

+0.06 

+0.49 

-62 

+0.10 

+0.010 

+0.013 

West 

APSF 

±0.07 

±0.06 

±45 

±0.40 

±0.005 

±0.004 


Amodel 

-0.07 

+1.66 

-117 

+ 1.23 

-0.003 

-0.003 


Asolar 

+0.16,+0.04 

+0.04,+1.07 

+85,-36 

-1.26,+1.55 

+0.005,+0.017 

+0.009,+0.018 


Aabun 

+0.18,+0.22 

+3.37,+0.37 

-63,+98 

+5.71,+2.14 

+0.016,-0.004 

+0.017, -0.002 


AnH 

±0.11 

±0.02 

±66 

±0.77 

±0.002 

±0 


ADistance 

±0.11 

±0.04 

±67 

±0.75 

±0.002 

±0.002 


AFI-BI 

-0.07 

+0.69 

-77 

+0.22 

+0.011 

+0.013 


(ANXB): vary the particle background component by 5%. 

(ACXB): fix the normalization of power law T = 1.4 at the expected value instead of letting it free to vary in the spectral fitting. 
(ACXB-T): fix the slope of power law at T = 1.5 and F = 1.3, respectively, for the CXB component. 

(ACXB-T (Extreme)): fix the slope of power law at F = 0.1 for the unresolved CXB component for the north pointing. 
(ASWCX): include solar wind charge exchange components in the model. 

(APSF): vary the mixing photons between each annulus by 5%. 

(Amodel): using a power law with three breaks (instead of two breaks) to model entropy profile. 

(Asolar): using solar abundance table Anders &; Grevesse (1989) and Lodders (2003) respectively. 

(Aabun): fix the metal abundance of the outermost bin at 0.1 Zq and 0.3^0 respectively. 

(AnH): vary the hydrogen column density by 20%. 

(ADistance): vary the redshift parameter by 5%. 

(AFI-BI): do not consider the difference in the responses between FI-CCD (XIS0 and XIS3) and BI-CCD (XIS1). 
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Fig. 13. — The Fe abundance profile of each direction. 
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Fig. 14.— Azimuthal scatters in temperature, density, entropy, 
and pressure calculated from eq. (10) for clusters studied in at 
least four directions with Suzaku. The results of other clusters are 
taken from the literature: Abell 1689 (Kawaharada et al. 2010), 
Abell 1246 (Sato et al. 2014), Abell 1835 (Ichikawa et al. 2013), 
PKS0745-191 (George et al. 2009), and the Perseus Cluster (Urban 
et al. 2014). Solid black cross: azimuthal scatter of simulated 
relaxed clusters; Dashed black cross: azimuthal scatter of simulated 
disturbed clusters (Vazza et al. 2011). 

dial profiles in the ICM properties are expected from 
such a system, any azimuthal spatial fluctuations present 
in the ICM reflect deviations from hydrostatic equilib¬ 
rium and/or a single-phase ICM. Numerical simulations 
demonstrate that the azimuthal variation of the ICM 
properties increases with radius beyond 0.5i?2oo (e.g., 
Burns et al. 2010; Roncarelli et al. 2013). For typi¬ 


cal cool-core clusters, which are known to have mostly 
regular X-ray images and appear to be approximately 
relaxed, recent X-ray observations reveal that approx¬ 
imate azimuthal asymmetry is commonly found in the 
outskirts of such clusters (e.g. Eckert et al. 2013; Walker 
et al. 2012a; Urban et al. 2014). 


6 .1. Azimuthal variations 

To quantitively describe the level of azimuthal varia¬ 
tion in RXJ1159+5531, we adopt the azimuthal scatter 
following Vazza et al. (2011): 


S'c(r) 


y- [Vi{r) ~ Y(r )} 2 
N Z+ [y (r) ]2 


( 10 ) 


where yi(r) is the radial profile of a given quantity in 
a section i. Y(r ) is the azimuthal average of this quan¬ 
tity, and TV is the number of azimuthal sections. Suzaku 
observed RXJ1159+5531 in four directions, TV = 4 and 
Y(r) — E^( r ) _ We obtained S C (R 2 00 ) of 0.076±0.050, 
0.086±0.060, 0.154±0.060, and 0.096+0.067 for the den¬ 
sity, temperature, entropy, and pressure, respectively, as 
shown in Figure [14} Its associated systematic uncertain¬ 
ties are listed in Thble 4. These measured values are com¬ 
parable to the predictions for relaxed systems (Vazza et 
al. 2011). Not many galaxy clusters have been observed 
with good azimuthal coverage near R vlI since the mod¬ 
erate spatial resolution of Suzaku restricts the observa¬ 
tions to nearby systems and their azimuthal studies can 
be observationally expensive. We found in the literature 
that four other galaxy clusters—-Abell 1689 (Kawaharada 
et al. 2010), Abell 1246 (Sato et al. 2014), Abell 1835 
(Ichikawa et al. 2013), and PKS0745-191 (George et al. 
2009 ^— also have been studied in four different direc¬ 
tions with Suzaku out to i? V ir- We compared the az¬ 
imuthal scatters in their gas properties at R 200 to that 
of RXJ1159+5531 in Figure [14] The nearby Perseus 
Cluster has been observed extensively with Suzaku in 
eight directions (still only a small fraction of its volume 
is covered). We also included the results of the Suzaku 
analysis of the Perseus Cluster (Urban et al. 2014). Az¬ 
imuthal scatters obtained with TV = 8 are expected to 
be ~20% larger than that obtained with TV = 4 (Vazza 
et al. 2011); thus we scaled the azimuthal scatters of the 
Perseus Cluster to its corresponding values for TV = 4 in 
Figure [14] The range of the azimuthal scatter of temper¬ 
atures is very large among these clusters. It is remark¬ 
able that the temperature scatter of RXJ1159+5531 is 
smaller than all the comparison clusters. In contrast, the 
range of the azimuthal density scatter is quite small and 
these clusters show comparable (and low level) density 
scatters (We note that the typical non-cool-core merging 
cluster Abell 1689, with the largest temperature scatter, 
has the smallest density scatter; thus the association of 
the density scatter with the cluster relaxation state is 
unclear). Consequently, the azimuthal variation in its 
entropy and pressure of RXJ1159+5531 at R 2 00 is the 


5 We note that the PKS0745-191 results given by George et al. 
(2009) may be compromised by inadequate background subtrac¬ 
tion, and, as argued by Walker et al. (2012a), its radial entropy 
profiles should start to deviate from the baseline at larger radii. 
However, it is unclear how this affects its azimuthal variation that 
we are interested in. 
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smallest among these systems. In Figure [14} we also over¬ 
plotted the predicted levels of azimuthaiscatters for re¬ 
laxed and disturbed clusters, respectively, based on simu¬ 
lations (Vazza et al. 2011). The lack of azimuthal scatter 
out to the R v i„ in the gas properties of RXJ1159+5531 
is a strong indication of it being a highly relaxed system 
(Roncarelli et al. 2013). 

In order to better understand the origin of the relaxed 
ICM, we inspected the Sloan Digital Sky Survey (SDSS) 
galaxy distribution on a scale of 20 Mpc within a \z\ < 
0.005 dispersion using spectroscopic redshift centered on 
RXJ1159+5531 as shown in Figure [l5| We marked the 
boundary of 3 Mpc centered on RXJ1159+5531. Some 
simulations show that galaxies start to be affected by ram 
pressure stripping within ~ 3 Mpc from the center of the 
cluster and the accretion flows of sub-structures also start 
around a few times R vlr (Cen et al. 2014). The surround¬ 
ing large-scale galaxy environment of RXJ1159+5531 in¬ 
dicates that the southwest directions seem connected to 
high density filamentary regions while the east direction 
is adjacent to “void” regions. However, there is no obvi¬ 
ous connection between its ICM properties around R v - lr 
and its large scale structures. 

Many other factors have been associated with large 
azimuthal scatter. These include unresolved substruc¬ 
tures (e.g., galaxies or subgroups, or even point sources 
for Suzaku), the presence of clumping gas or relativis¬ 
tic plasma, and gas displacement associated with bulk 
motions. The strikingly symmetrical gas properties 
found in RXJ1159+5531 place constraints on the im¬ 
portance of the above processes. The lack of significant 
azimuthal scatter in the gas properties out to R v ; r of 
RXJ1159+5531 is strong evidence of it being a highly 
evolved system and hydrostatic equilibrium to be a very 
good approximation. 

6 .2. The origin of the unexpected entropy profile 

As discussed in the Introduction, cluster entropy (K = 
kT/n^ 3 ) profiles at large radii have been frequently ob¬ 
served to deviate from the power-law profile (K oc r 11 ) 
predicted by “adiabatic” cosmological simulations that 
account only for gravity (Voit et al. 2005), and several 
explanations have been proposed to reconcile this dis¬ 
crepancy. 

The proposed mechanisms all depend on the dynamical 
state or mass scale of the cluster, which provides us an 
opportunity to place observational constraints on them 
by comparing the entropy profiles of galaxy clusters of 
different masses and at different evolutionary stages. For 
example, gas dumpiness has been predicted to be more 
prominent among more massive galaxy clusters. This 
is because more massive clusters formed more recently 
and are more perturbed (Vazza et al. 2013; Battaglia 
et al. 2010), although Nagai & Lau (2011) argue this 
is because lower mass systems have a larger fraction of 
lower-temperature gas, and cooling removes high-density 
gas out of the detectable energy range. The presence 
of non-thermal pressure support (associated with depar¬ 
tures from hydrostatic equilibrium) is also expected to be 
more pronounced in dynamically younger (and usually 
more massive) clusters (Shi & Komatsu 2014). Avestruz 
et al. (2014) show that the non-thermal equilibrium be¬ 
tween ions and electrons is more significant in hotter 
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Fig. 15.— Galaxy number density map coded with their r-band 
luminosities from SDSS catalogue centered on RXJ1159+5531. 
The spectroscopic redshifts of these galaxies are used for selection 
and only those within a 1500 kms -1 dispersion are included. The 
white circle indicates the physical radius of 3 Mpc. 

clusters due to their longer equilibration time. Fujita 
et al. (2013) also found that cosmic ray acceleration is 
more effective in hotter clusters. On the contrary, the 
scenario proposed by Cavaliere et al. (2011) predicts 
that the entropy drop should be minimal in dynami¬ 
cally young systems since merger shocks are more active 
among them and highly evolved clusters should deviate 
more from the adiabatic ~ r 11 profile than most other 
clusters. The entropy at R 200 measured in all directions 
of RXJ1159+5531 are consistent with the theoretical ex¬ 
pectations for purely gravitational processes operating 
in galaxy clusters. Its enclosed gas fractions within R v j r 
derived for all directions are consistent with the cosmic 
value. Compared to other massive clusters observed out 
to i?vir by Suzaku , RXJ1159+5531 is less massive and 
highly evolved. The ICM properties of RXJ1159+5531 
thus disfavor the “weakening of accretion shock” expla¬ 
nation (Cavaliere et al. 2011) and are consistent with 
other mechanisms, though they do not require them ei¬ 
ther. 

Deeply related to these mechanisms, the entropy be¬ 
havior in the cluster outskirts has been frequently as¬ 
sociated with the large scale structure. For example, 
using SDSS, Kawaharada et al. (2010) found that in 
Abell 1689 the direction adjacent to higher galaxy den¬ 
sity on the large scale shows larger entropy than in other 
directions adjacent to smaller galaxy density regions. 
They propose that the thermalization of the ICM oc¬ 
curs faster along overdense filamentary structures than 
along low-density void regions. However, Kawaharada 
et al. (2010) employed photometric redshifts when they 
construct the large scale galaxy density map; the un¬ 
certainty of photometric redshifts can be very big. Ur¬ 
ban et al. (2014) studied the Perseus Cluster in eight di¬ 
rections with extensive Suzaku observations. They also 
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found that directions along the major axis show higher 
entropy than other directions. They propose that more 
frequent major mergers occur along the major directions 
and such mergers can effectively destroy clumping gas, 
although we note those directions with high entropy at 
i? 2 oo °f the Perseus Cluster are associated with high 
temperatures instead of low density gas. However, we 
measure the same entropy in each direction at i? v ir for 
RXJ1159+5531, while its large scale galaxy distribution 
is asymmetric (Figure |l5|). Such a lack of connections be¬ 
tween gas properties in the outskirt of a cluster and the 
large scale structures has also been reported for Hydra 
A (Sato et al. 2012) and Abell 1246 (Sato et al. 2014). 

Finally, it has been suggested that in evolved clusters, 
the low-density magnetized ICM will develop a Magne- 
tothermal Instability (MTI), especially in the outer re¬ 
gions (e.g., Parrish et al. 2008). The MTI drives tur¬ 
bulent motions and deviations from hydrostatic equilib¬ 
rium in the cluster outskirts. Since our measurements 
for RXJ1159+5531 indicate that it must be very close to 
hydrostatic equilibrium, we can infer that the MTI must 
not be operating effectively in this system. 

7. SUMMARY 

We revisited the Chandra and Suzaku observations of 
the fossil-group/poor-cluster RXJ1159+5531. We pre¬ 
sented the analysis of the Suzaku observations of this 
system of three new directions. We also revised the pre¬ 
vious study with deep Chandra ACIS-I observations cov¬ 
ering out to the virial radii of the north direction of this 
system. Below, we list a few highlights of this paper: 

• We have resolved more than half of the CXB compo¬ 
nent into point sources with deep Chandra exposure in 
the north direction. With the refinement of the Chandra 
data, we measured an entropy profile rise all the way to 
R VII and a baryon fraction within Rios consistent with 
the cosmic value, confirmed our previous study of the 
this system in the north direction. 

• The gas properties (temperature, density, entropy, 
and entropy) of all four directions (north, south, east, 
and west) have very similar values. The azimuthal scat¬ 
ter in entropy is much smaller than that found in other 
galaxy clusters. The good agreement of the entropy with 
that predicted from gravity only simulations at R 2001 
coupled with the small azimuthal variations found for 
the ICM spectral properties, especially for the entropy, 
indicates that the ICM of RXJ1159 is very close to hydro¬ 
static equilibrium and, as such, unusually relaxed com¬ 
pared to other clusters. 

• Since the virial entropy in all directions are consis¬ 
tent with the theoretical expectation within uncertain¬ 
ties and the gas fractions derived for all directions are 
consistent with the cosmic value, there is no need to 
invoke gas dumpiness or any non-gravitational process 
in RXJ1159+5531. This result is inconsistent with the 
weakening of accretion shock explanation given by Cav- 
aliere et al. (2011) although their model fits the entropy 
profiles of most cool-core clusters (Walker et al. 2013). 

• We inspected the galaxy density distribution cen¬ 
tered on RXJ1159+5531 on a scale up to 20 Mpc. The 
gas properties near R v - lr seem unrelated to its large struc¬ 
ture environment. 
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TABLE 4 

Azimuthal scatter of gas properties at R200 with their systematic budgets 



Temperature 

Density 

Entropy 

Pressure 

best-fit 

0.076 ±0.050 

0.086 ± 0.060 

0.154 ±0.067 

0.096 ±0.067 

ANXB 

±0.019 

±0.038 

±0.052 

±0.057 

ACXB 

±0.062 

±0.015 

±0.048 

±0.063 

ACXB-r 

±0.106 

±0.035 

±0.147 

±0.074 

ASWCX 

±0.025 

±0.020 

±0.028 

±0.021 

APSF 

±0.094 

±0.033 

±0.154 

±0.088 

Amodel 

±0.024 

±0.022 

±0.042 

±0.007 

Asolar 

±0.061, ±0.013 

±0.040, ±0.024 

±0.073, ±0.044 

±0.028, ±0.050 

Aabun 

±0.066, ±0.024 

±0.058, ±0.018 

±0.073, ±0.020 

±0.072, ±0.077 

AnH 

±0.061 

±0.031 

±0.050 

±0.105 

ADistance 

±0.052 

±0.007 

±0.015 

±0.021 

AFI-BI 

±0.047 

±0.018 

±0.026 

±0.063 

Relaxed 

0.091 ± 0.097 

0.081 ± 0.061 



Perturbed 

0.166 ±0.277 

0.128 ±0.071 




Systematic error definitions are the same as for Table 3. The la statistical error on the best fit is obtained from the standard deviation 
of measurements from 10000 sets of Monte Carlo simulations of the best-fitting spectral models. Predictions from numerical simulation of 
azimuthal scatters in temperature and density at R200 i n relaxed and perturbed clusters are also listed (Vazza et al. 2011). 
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